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We evaluate, analytically, a specific class of eighth-order and tenth-order QED contributions to 
the anomalous magnetic moment of the muon. They are generated by Feynman diagrams involving 
lowest order vacuum polarization insertions of leptons I = e,fi, and r. The results are given in the 
form of analytic expansions in terms of the mass ratios m e /m M and m M /m T . We compute as many 
QO ' terms as required by the error induced by the present experimental uncertainty on the lepton masses. 

We show how the Mellin-Barnes integral representation of Feynman parametric integrals allows for 
an easy analytic evaluation of as many terms as wanted in these expansions and how its underlying 
algebraic structure generalizes the standard renormalization group properties. We also discuss the 
» I , generalization of this technique to the case where two independent mass ratios appear. Comparison 

C*l |i with previous numerical and analytic evaluations made in the literature, whenever pertinent, are 

also made. 



PACS numbers: 



INTRODUCTION 



The present experimental world average of the anomalous magnetic moment of the muon a^, assuming CPT- 
variance, viz. a M + = a M - , is 

CM 

> 1 a(° xp) = 116 592 080 (63) x 10~ n (0.54 ppm) , (1.1) 

oo ; 

where the total uncertainty includes a 0.46 ppm statistical uncertainty and a 0.28 ppm systematic uncertainty, 
combined in quadrature. This result is largely dominated by a series of precise measurements carried out at the 
Brookhaven National Laboratory (BNL) during the last few years, by the E821 collaboration, with results reported in 
ref. [l|| and references therein. The prediction of the Standard Model, as a result of contributions from many physicists 

,(SM) = i 1fi KQ1 78K ffi-n v 10- 11 



aj° M > = 116 591 785 (61) x 10" 11 , (1.2) 

where the error here is dominated at present by the lowest order hadronic vacuum polarization contribution uncertainty 
(±46.6 x 10 -11 ), as well as by the theoretical uncertainty in the hadronic light-by-light scattering contribution, 
estimated to be ±40 x 10 -11 . Errors here have also been combined in quadrature. The results quoted in Q1.1J1 and 
(|1.2j) imply a 3.4 standard deviation in the difference 

4° XP) " a » M) = ( 295 ± 88 ) x 10_n ■ ( L3 ) 

This 3.4 a deviation deserves attention. Ideally, and before one can attribute the present discrepancy to new 
Physics, one would like to reduce the theoretical uncertainties as much as possible, parallel to a new experimental 
effort toward an even more precise measurement of a M 0, E|. It is also important to reexamine critically the various 
theoretical contributions to Eq. l|1.2j) : primarily the hadronic contributions of course, but also the higher order QED 
and electroweak contributions. It would be reassuring to have, at least, two independent calculations of some of these 
contributions, as well as of the higher order estimates. The purpose of this article is a first step (albeit a small one) 
in this direction. 
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We shall be concerned with a specific class of eighth-order and tenth-order QED contributions to illustrated by 
the Feynman diagrams shown in Fig. 1 and Fig. 2 below. The contribution from the diagrams in these figures with 
no r vacuum polarization insertions i.e., diagrams (A), (B) and (C) in Fig. 1 and diagrams (A), (B), (C) and (D) in 
Fig. 2, have been estimated numerically by Kinoshita and Nio in refs. [5|| and [6j (see also ref. [?(]). There are also 
analytic results of the same diagrams, published by Laporta [1, [3], in the form of analytic expansions in terms of the 
electron to muon mass ratio To our knowledge, the diagrams in Fig. 1 (D) and Fig. 2 (E) involving a tau loop 
vacuum polarization insertion have not been considered before. We find for these contributions the following results: 



where the errors are the ones induced by the present accuracy in the determination of the lepton masses. As one can 
see from the other results reported in the Appendix, the contribution from the diagrams in Fig. 2 (E) is of the same 
size as the one from the diagrams in Fig. 2 (D). 

In a previous paper [13], it has been shown how the Mellin-Barnes integral representation of Feynman parametric 
integrals combined with the converse mapping theorem [ll[ , allow for an easy evaluation of as many terms as wanted 
in the asymptotic behaviour of Feynman diagrams in terms of a mass ratio. We shall show how this technique applies 
to the evaluation of the class of diagrams considered here, and how to generalize it to the case where two independent 
mass ratios, like ^ and ^f- in the diagrams in Fig. 1 (D) and Fig. 2 (E), enter into consideration. We shall do that 
in a rather detailed way for various reasons: 

• We find, regretfully, that many analytic calculations found in the literature provide very few details about the 
techniques employed; yet it would be useful to compare the efficiency of different methods when planning new 
more complex calculations. 

• Giving details on some of the intermediate steps of complex calculations can only help to see the qualities and 
limitations of the methods employed and, occasionally, to spot potential errors. 

• As explained in ref. [13] many years ago, there is a large class of higher order contributions which can be estimated 
using renormalization group arguments: powers of log terms at a given order are algebraically related to lower 
order contributions. Subsequent applications have been made by many other authors [32]; however, very little is 
known about how to extend renormalization group arguments to subleading terms: constant terms and powers 
of log 2a suppressed by ^f- powers. The Mellin-Barnes technique that we are advocating provides a precise 
answer to that question: as we shall see, in the Mellin-Barnes representation, the asymptotic contributions at 
a given order factorize in terms of well defined moments of lower order contributions. 

• Having a powerful technique to obtain asymptotic expansions provides a useful alternative to the computation 
of exact analytic expressions, which are often very complicated and cumbersome. After all, ratios of masses are 
known from experiment only to a fixed accuracy; the associated error propagates into a numerical uncertainty of 
the exact analytic result in any case. Computing as many terms in the corresponding asymptotic expansion, as 
required by the experimental precision in the masses involved, provides an easier alternative to the computation 
of an exact analytic result, with the same practical accuracy. The calculations in this paper have all been done 
within this spirit. 

• Finally, we think that it is important to have an independent way to check the precision of the contributions 
evaluated numerically. Many of the multidimensional integrals involved in higher order calculations are far from 
trivial, which has obliged the experts to develop skillful methods. The numerical results are often dominated by 
a statistical error which is larger than the error induced by the experimental determination of the lepton masses. 
As we shall see, within this limited precision, we find that all the results of Kinoshita and collaborators [3, [3, 0] 
checked in this paper are correct, within less than one standard deviation of their estimated error, which is a 
remarkable performance. 

The paper has been organized as follows. The next section summarizes well known results about lepton vacuum 
polarization and the way we choose to express the various contributions to the muon anomalous magnetic moment. 
Section III gives a detailed description of the Mellin-Barnes technique and the converse mapping theorem as applied 
to the eighth-order and tenth-order contributions that we are considering. We also discuss in this section the 
underlying algebraic properties which generalize the usual renormalization group constraints previously discussed in 
the literature. Section IV contains a discussion of the evaluation of the various moment integrals which appear in 
the intermediate steps. They can all be done analytically and we give explicit expression for all of them; but we also 
show how in practice, the evaluation of a finite number of terms in the ultimate expansion, only requires the explicit 




and 




(1.4) 
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knowledge of approximate expressions for these moments. This is important to know, in view of more complicated 
situations that one may very well encounter when considering other Feynman diagrams. The final evaluation of the 
eighth-order diagrams is done in Section V and of the tenth-order diagrams in Section VI. Section VII is dedicated 
to the generalization of the techniques discussed in previous sections to the case where three mass scales appear. The 
calculation of the diagrams in Figs. 1 (D) and 2 (E) are non-trivial examples of applications of this new technology. 
Finally, we collect in an Appendix all the results that we have computed in this paper. 



II. PHOTON PROPAGATOR AND THE MUON ANOMALY 



We are interested in the replacement of a free photon-propagator by the dressed propagator 

q a qp\ i i 



D a /3(q) 



i 9a(5 



9 2 i + E/n('V) 



(2.1) 



where 11^ (q 2 ) denotes the proper vacuum polarization self-energy contribution induced by a lepton loop I = e,fj,,r 
and a is a parameter reflecting the gauge freedom in the free-field propagator (a = 1 in the Feynman gauge). In fact, 
all the diagrams which we shall be considering here are gauge independent and, therefore, the term i(l — a)^|^- does 
not contribute to their evaluation. 

The perturbation theory expansion generates a series in powers of the functions 11^ (q 2 ). Here we shall be concerned 
with a particular selection of the terms in this series, namely those with three and four powers of the lowest order 
n( z )(q 2 ) self-energies; more precisely with the terms: 



l + Ejiicv) 



- 3 



n ( e ) (g 2 )l (q 2 ) _ 3 n( e ) (q 2 ) few ( q 



+ 



u^( q 2 )} nM(« 2 ) 



n (e) (q 2 )l rr>V) + 6 [n (e V)l [n ( ^V) 



4 n (e V) 



ir>V 



+4 



IT^V)] IT>V) 



(2.2) 

(2.3) 
(2.4) 
(2.5) 



When inserted in the free-photon propagator of the lowest order one loop muon vertex, these terms generate the 
four types of Feynman diagrams collected in Fig. 1 and the five types of diagrams collected in Fig. 2. The Feynman 
diagrams in these figures give, respectively, eighth order and tenth order contributions to the anomalous magnetic 
moment of the muon, which are enhanced by powers of log ^f- factors. This is why we are interested in these terms. 

The reason why we only keep the terms in (|2.3p and (|2.5p with one power of H^(q 2 ) is because, in spite of the 

enhancement by log ^f- factors, the r-loop induces a suppression factor ~ ^f. Contributions from two or more 

powers of Il( r )(g 2 ) factors to the muon anomaly will be, therefore, neglected. 

A very convenient integral representation for the contribution to the muon anomaly from the graphs in Fig. 1 and 
Fig. 2 can be obtained as follows. First we write a dispersion relation for the on-shell renormalized photon propagator 
induced by electron self-energy loops only. This generates the following set of relations: 

dt q 2 ( Am 2 \ 

Pi — - , J = 1,2,3,4, (2.6) 



n( e V 



t t-f 



with 



Pi 



P-2 



P3 



Pi 



4m 2 \ 


1 


. ) 


7T 


4m 2 \ 


1 


. J 


7T 


4m 2 \ 


1 


. J 


7T 


4m 2 \ 


1 




7T 



-i m n (e) (t), 

- {2 ReIT (e) (t) Imn (e) (t)} , 

Ren (e) (t)l 2 ImIT (e) (t) - [lmIT (e) (t) 
Rell^ (t)l 3 Imn< e ) (t) - 4Ren( £ ) (t) [imll^ (t) 



(2.7) 
(2.8) 
(2.9) 
(2.10) 
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Using these representations, the contribution to the muon anomaly from the various terms in the series in 
Eqs. I|2.2p . I|2.3p and (|2.4p . (|2.5p can be viewed as the convolution of electron spectral functions, modulated by 
powers of the muon or tau self-energy functions, with a free-photon propagator replaced by a fictitious massive 
photon propagator: 



1 

T 



(-1) 



t + ie 



(2.11) 



The muon vertex loop integral over the virtual g-momenta can then be traded by a Feynman rr-parameter integral 
(see the review article in ref. [l5| and the earlier references therein), with a net contribution to the muon anomaly, 
up to an overall combinatorial factor Fr^ p \ which can be read from the expansion in Eqs. (|2.2p to f|2 . 5[) : 



Aj,p) _ H / * 



dx- 



x 2 (\~x) 



Ft 



1 "V 

1 — x p 



pA — 



(2.12) 



where p = 0, 1, 2, 3 when I — /i and p = 1 when I — r, while the index j counts the number of electron loops in the 
vacuum polarization; therefore 



Ft 



(3,0) 



1, F ( 2,i)=3, F (1 , 2) =3, 



(2.13) 



and 



F 



(4,0) 



1 



F(3,l 



(2,2) 



-6, F (1 



(1.3) 



(2.14) 



Notice that in the process of trading the integral over the virtual g-momenta by the Feynman parameter x, one has 
obtained the effective replacement: 



(2.15) 



TO,, . 

1-x 



in the muon and tau self-energy functions. 

The lowest order vacuum polarization self-energy functions in QED are well known. We shall use the representations 
given in ref. [l6j]. With 



5 = \ 

the lowest order spectral function for the electron is 



4to2 



-Imn>) (t) = -s(-- -S 2 ) 9 (t - 4m; 



(2.16) 



(2.17) 



and the real part [5 = Jl — 



Am'l 



IS 



Ren 



(e) 



8 Iri! r/ 1 IraA, I 1 " 5 ! 



(2.18) 



We shall also need the expression for the vacuum polarization self-energy induced by a muon loop, in the euclidean 
region, and as a function of the Feynman x-parameter: 



n 



(/<) 



to,, 

1-x " 



5 4 4 
9 + 3a? ~ 3X 2 " 



4 \ 



3a; 3 J 



log(l - x) 



(2.19) 



as well as the one induced by a tau loop which, for our purposes, it is more convenient to keep in the Feynman 
parametric integral representation: 



n(r) 



TO,, 

1-X 



dz2z(l — z) log 



1-x 



rat 



(2.20) 
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III. THE MELLIN BARNES REPRESENTATION AND THE RENORMALIZATION GROUP 



It is now useful to introduce the two Mellin-Barnes integral representations [3c 

C s +20O 



X 2 (l 



1 



ds 



Ami 



4to; 



x 2 + -\(l-a 

* c s —ioo 

where the threshold scale Am 2 is chosen for later convenience, and 134 



x 2s (i- x y- s r(s)r(i- s ) . 



(3.1) 



log 



1 



1 



tot, 

-z(l-z)-^ 

x mi 



c±-\-ioo 




-z(l-z) 



™r(i- t ). 



where the integration paths along the imaginary axis are defined in the fundamental strips flll | : 

c s = Re(s) € ]0, 1[ and c t = Re(t) £ ] - 1,0[. 

The interest of these representations lies in the property that the dependence on the physical mass ratios 



(3.2) 

(3.3) 
L and 



— f is then fully factorized and the Feynman parametric integrals one is left with are then those of a massless case. 
The only new feature is that they have to be computed as functions of the Mellin s-complex variable, in the case of 
electron and muon loops only; and as functions of the (s , i)-complex manifold in the case of electron, muon and tau 
loops. In all the cases we shall be concerned with here, it is possible to obtain analytic expressions for these functions 
in a rather straightforward way. In fact, they turn out to be rational functions of products of Gamma functions and 
Polygamma functions which depend linearly on the Mellin variables s, or s and t. 

Let us first discuss the case where we only have one Mellin s-complex variable. It corresponds to the Feynman 
diagrams where only two mass scales are involved (i.e. the case where I = fi in Eq. Ij2.12p ). Using the representation 
in Eq. (|3.1|l one can rewrite Eq. I|2.12p . for a fixed p and j, as follows 



a 1 
7T 2?ti 



c s +ioo 



ds 



Ami 



M(s) 



(3.4) 



with 



M(s) = T(s)T(l - s) / dxx 2s (l-x) 



F, 



U,p) 



-TO, 



00 dx 
t 



Amt 



Ami 



(3.5) 



The converse mapping theorem relates the asymptotic behaviour of a M as a function of the small mass ratio Am^/m 2 ^, to 
the singularities of the integrand M(s) as a function of the Mellin s-complex variable. For m 2 <§; to^ the appropriate 
s-singularities are those in the left-hand-side of the fundamental strip and they are all on the negative real axis. The 
precise relation goes as follows: with p € R and k GN, the function Ai(s) in the l.h.s. of the fundamental strip has a 
singular expansion of the type (ordered in increasing values of p) : 



(S + P ) fe +! 



(3.6) 



The corresponding asymptotic behaviour of (ordered in increasing powers of p) is then: 



a^ v (-l) fc 



4m„ 



log' 



4to„ 



(3.7) 



and the problem is then reduced to the calculation of the a Pi fc residues. This, in turn, requires the evaluation of the 
moment integrals 



00 dx 
t 



Ami 



Pi 



( Ami 



V t 



and 



dx x^il-x) 1 



(3. 
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which appear in Eq. Q3.5P 

Let us consider a little bit further these two types of moments. First, the moments of the spectral functions 
Pj (nr*) defined in E Q S - EH to IpTO)) i.e., 



(;)'*<•> 



t V t 



ft 



4mj 



(3.9) 



One can immediately see that these moments can also be seen as the Mellin transform of the spectral functions Pj 



with respect to the variable £ = since 



4m„ 



(3.10) 



where in fact the integral over £ only runs from zero to one, because of the threshold factor 9{t—Am 2 e ) in the lowest order 
spectral function in Eq. l|2.17p . The interest of this point of view lies in the fact that the Mellin-Barnes representation 
relates the singular expansion of any of these Rj(s) moments to the asymptotic expansion of the spectral functions 
Pj(^). Indeed, from Eq. I|3.10p and the asymptotic behaviour of pj(£) at small £ and at large £, there follows that 



r s +200 



Pj 



(€)= 0'i / dS Q) R ^ S)l ^ r «= Re (s) G ]0,oo[. 



(3.11) 



In other words, the direct mapping (in the sense of ref. [ll|) establishes a precise relation between the asymptotic 
expansion of the spectral functions pj(£) at small £ (which is known from their explicit analytic expression) and the 
singular series of the moments Rj(s), for s < 0, that we are interested in. We shall come back later on to these 
relations in somewhat more detail. 

The second type of moments in Eq. I|3.8p are also well defined Mellin transforms; in this case with respect to the 
invariant photon momenta 



1-x 



(3.12) 



flowing in the basic muon vertex diagram. Indeed, one can easily verify that 



/ Ot\P 

(-) "pW 



dx x 2s {\ - x) 1 



du> lu s 1 



ir» f ^- x ml 



4 + W \ y/4 + UJ + 



nW(-i 



(3.13) 



Again, the evaluation of the singular series associated to the moments fi p (s), for s < 0, is very much facilitated by 
the direct mapping which relates them to the asymptotic expansion of lj m 2 l )\ P for small w. 

The remarkable property of the Mellin-Barnes representation is precisely the factorization in terms of moment 
integrals as shown in Eq. I|3.5p . It is in fact this factorization which is at the basis of the renormalization group 
properties discussed by many authors. The algebraic factorization above, however, is more general and it also shows 
the full underlying renormalization group structure at work. The classical renormalization group constrains which have 
been exploited in the literature only apply to the evaluation of leading asymptotic behaviours {powers of logarithms 
and constant terms). In the example above, this is encoded by the properties of the Mellin singularity at s — 0. This 
singularity governs the spectral function moments Rj(0) on the one hand, as well as the contribution to the muon 
anomaly from vacuum polarization muon loop insertions (p loops): the moments f2 p (0). Notice that the moments 
Rj(0) are UV-singular, hence the need of the s- regularization. The singularity is related to charge renormalization 
and therefore to the QED /3-function at a given order in perturbation theory [12j, in our case the lowest order. The 
predictive power of the renormalization group structure lies in the fact that once we know the two types of moments 
in Eq. (|3.8p at a given order in perturbation theory we have a prediction at a higher order for the convolution which 
in our case gives a M in Eq. I|3.4p . What is new here is that this factorization extends as well to the subleading terms 



in the expansion modulated by inverse m^/mg-powers. What governs these terms is nothing but the residues of the 
successive Mellin singularities in the negative real axis. Here, the expansion in inverse m^/m^-powers is analogous 
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to the l/Q 2 -expansion in the Operator Product Expansion of Green's functions in QCD, while the moments Rj(s), 
for s < 0, are the equivalent of the so called vacuum condensates. The big difference, of course, is that in QED these 
condensates are moments of spectral functions which are known at all values of t and, therefore, can be calculated 
explicitly at a given order in perturbation theory. Notice, however, that also in QED these moments (condensates) 
are a priori ill defined because they are singular. The Mellin-Barnes technique provides a systematic regularization 
and hence a precise separation of short-distance and long-distance effects. 

The case of two Mellin (s , i)-complex variables, corresponding to Feynman diagrams with both electron loops and 
tau loops, is conceptually more complicated and it requires a specific treatment which we shall provide in Section VII. 
A detailed discussion of the very interesting underlying mathematics which governs this case will be given in a 
forthcoming publication (2lj . 

We finally give below the exact explicit expressions, in the Mellin-Barnes integral representation, of the contribution 
to the muon anomaly from each of the specific set of Feynman diagrams in Fig. 1 and Fig. 2. 

A. Eighth Order Contributions 



X 




(D) 



Fig. 1 Eighth-order Feynman diagrams with lowest order vacuum polarization electron-loops and a r-loop which contribute to 
the Muon Anomaly and are enhanced by powers of log — a factors. The dots indicate the other diagrams with different 

permutations of the lepton-loops. 

• Three Electron Loops, Fig. 1(A) [one diagram]: 



^ ee) = © 4 i / K^f) St ^--s)^(s)R 3 (s). (3.14) 

c s — too 
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Two Electron Loops and One Muon Loop, Fig. 1(B) [three diagrams]: 

c s +ioo _ s 

4 eeAl) = (-) 4 A / ds(^4) r( s )r(i- s )^( s )i? 2 ( s ). (3.15) 

p \7T/ 27TI J \ m~ J 

c s —ioo 

One Electron Loop and Two Muon Loops, Fig. 1(C) [three diagrams]: 



C s +'iOO _ 

a\ 4 3 f , /4m 2x 



4 ew0 = (r ^- / * -r r( s )r(i- s )n 2 ( s )i?x( s ). (3.16) 



• Two Electron Loops and One Tau Loop, Fig. 1(D) [three diagrams]: 

Using the representation in Eq. (|2.20p for the tau self-energy function as well as the Mellin-Barnes representation 
in Eq. 113.21) one easily gets the expression 



c 3 +ioo c t +ico , x -t 

ioo 

r(t) n 



Ct — 200 



r( s )r(i - fl ) -^-r(i - t) e(«, t) J? 2 (s) , (3.17) 

where 0(s,i) is the Feynman parametric integral (see Eqs. (|3.ip and l|3.2p ) 

0(«,t) = / dxx^il-x) 1 - 8 [ dz (-^—-z{l - z) 



Jo 



_ r(2 - t)r(2 - 1) r(i + 2 S - 2t)r(2 - g + 1) 

- ( ~ 2) r(4-2t) r(3 + «-t) ■ (3 ' 18) 

Notice that the dependence on the variables s and t in the function 0(s,t) is not factorized. It is this fact that 
requires new technical considerations which we shall discuss in Section VII. 

B. Tenth Order Contributions 

• Four Electron Loops, Fig. 2(A) [one diagram]: 



a (eeee) 



c s +zoo 

,2- 



Gr) 8 ^ / dS (^) r(8)r(1 - s)fio(8)J24(s) ' (319) 

c s — zoo 

Three Electron Loops and One Muon Loop, Fig. 2(B) [four diagrams]: 



(eee/j) 

It 



c s +ioo 

.2 



(-) 5 (-4) 4t^ / r( S )r(l- S )Ox( S )i? 3 ( S ). (3.20) 

\7T/ Z7TI J \ / 
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x 




Fig. 2 Tenth-order Feynman diagrams with lowest order vacuum polarization electron-loops and a r-loop which contribute to 
the Muon Anomaly and are enhanced by powers of log ^f- factors. The dots indicate the other diagrams with different 

permutations of the lepton-loops. 

• Two Electron Loops and Two Muon Loops, Fig. 2(C) [six diagrams]: 



c s +ioc _ 

^^Gf)" fcrf / KSjf) S nmi-s)n 2 (s)R 2 (s). (3.21) 

C a — lOO 

• One Electron Loop and Three Muon Loops, Fig 2.(D) [four diagrams]: 



c s +ioc _ 

a (e^) = ^ 5 (_^) J ds (^j T{s)T{l-s)^{s)R 1 {s). (3.22) 

c s — too ^ 

• Three Electron Loops and One Tau Loop, Fig. 2(E) [four diagrams]: 

The only change here, with respect to the representation in Eq. (|3. 1 7|1 . is the convolution with the moment 
Rz(s) (instead of Rz{s)) and the combinatorial factor (—4) [instead of (3)]: 
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c 3 +ioo e t +ioo , \ — t 



r( s )r(i - s) -^r(i - t) e(«, t) r 3 (s) , (3.23) 



where Q(s,t) is the function defined in Eq. (|3- 18|1 . 



IV. THE MOMENT INTEGRALS 

We are now left with two types of moment integrals: the Rj(s) moments of the electron-positron spectral functions 
in Eq. (|3.10p . (the integrals over the variable £, or equivalently the variable 6 in Eq. (|2.16p ): and the O p (s) moments 
of the euclidean muon vacuum polarization, (the integrals over the invariant photon momenta u, or equivalently 
the Feynman x-parameter) . It is a luxury that these integrals can all be done analytically. Here we shall give the 
results and comment on a few technical points; in particular on their relation to the underlying renormalization group 
properties. 

• The O (s) Moment. 

This moment, which appears in a^ eee ^ and a^ eee \ corresponds to the trivial x-integral 

^ / x f 1 , , J x 2 \ s r(i + 2s)r(2 - s) , . 

The particular value f2o(0) = \ is precisely the coefficient of the lowest order muon anomaly, the Schwinger 
term: a M = (^) [Oq(0) = 41. Notice that flo(s) is IR-singular at s = — 1/2, — 3/2, ■ • • . It is precisely these 
IR-singularities that are at the origin of the contributions in odd powers of m e /m^ in a^ e \ a^ ee ^ and 
^(eeee) one encounters in the explicit calculations. 

• The Rx(s) Moment. 

This moment appears in the expressions for ajf^ and in Eqs. (|3.16p and (|3.22|) and it corresponds to 



the integral 



R^s) = J dS 26(1 - 5 2 r 1 \& (l - l/^j , (4.2) 



which can be trivially computed, with the simple result: 

i r(2 + s ) 

explicitly showing the singularities in the negative real axis (the origin s — included). Notice that, because of 

r / 2\ p- 

the zeros provided by the factor jrpq~^Y) there cannot be terms of O [^f J for p = 5/2 + n, n = 0, 1, 2, • • • 
in vacuum polarization contributions to a M from one electron loop. 

The two moments flo(s) and Ri(s) are the only quantities required to evaluate the contribution to affi from the 
fourth order Feynman diagram in Fig. 3, which, in our representation, is given by the integral 

< } = Q 2 oh J ds {^) * nmi - *) n (s) r,(s) ■ (4.4) 

c 8 — ioo 

Since the moments n (s) and Ri(s) are known analytically, the singular series expansion of the integrand can be 
computed up to as many terms as one wishes; e.g., 
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X 




e 



Fig. 3 Fourth-order Feynman diagram with a lowest order vacuum polarization electron-loop. 




-2(TT7F + (i- log2 



') 



l + s 



1 



+ 



(4.5) 



and the converse mapping theorem tells us how to read in a straightforward way the corresponding asymptotic 
expansion contribution [35]: 



The reason why we consider this well known contribution here is that it allows us to discuss the factorization properties 
of the Mellin-Barnes representation and its relation to the underlying renormalization group structure in a case which, 
in spite of its simplicity, illustrates the generic features rather well. Notice that the singular expansion in Eq. (|4.5p 
results from the combination of the Laurent series of three factors: the geometric series factor r(s)T(l — s), which 
appears in the original Mellin-Barnes representation in Eq. (|3,lj) . and the two moments flo(s) and Ri(s). The residue 

of the leading singularity in 1/s 2 , which provides the coefficient of the log term in Eq. (|4.6j) . can be read off directly 

from the asymptotic behaviour of the lowest order spectral function p {^—^j ~ §, (which in turn is correlated to 

the residue of the leading term in the singular expansion of Ri(s) x |~) and the lowest order result Qq = 1/2 ; 
the geometric factor T(s)r(l — s) providing the extra 1/s factor. This is precisely the leading prediction of the 
renormalization group in this case [12 ] . 

The next -to-next-to-leading term in the asymptotic expansion in Eq. (|4.6p (the second line) is governed by the 
s+1 / 2 term in the singular expansion in Eq. I|4.5p (the second line in this equation). The origin of this singularity 
is the moment ^o(s), which is singular at s = —1/2; the other factors are regular: T(— 1/2)T(1 + 1/2) = —ir and 
i?i(l/2) = —j. Again, the residue of the 3+ \/ 2 -singularity of CIq(s) can be read off from the leading term in the 
asymptotic expansion of the integrand function in Eq. (|3 . 13[) i.e., 




(4.6) 




(4.7) 
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We see, therefore, how the factorization of the Mellin-Barnes representation allows one to fix a non-trivial coefficient 



of the asymptotic expansion of an O 



quantity from the knowledge of two numbers which appear at a lower 



O (^) and have been very easily identified. 

What about the next-to-leading term? This refers to the factor — 1| in the first line of Eq. l|4.6p , which is not fixed 
by simple renormalization group arguments. The reason why simple renormalization group arguments fail to fix this 
term is due to the fact that it does not originate from products of leading terms of singular expansions only. More 
precisely, the coefficient of the \ term in Eq. (|4.5[) also depends on the next -to-leading terms of the Laurent series 
of each of the individual factors fla(s) , and Ri(s), These terms, however, can also be easily obtained without having 
to calculate explicitly the functions Oo(s) and R\(s). Indeed, since fio(s) is regular at s — ► we can simply Taylor 
expand the integrand and find 



n ( s ) 



1 

'o 2 



■0(s 2 



(4i 



By contrast, the moment Ri(s) is singular at s — > 0; but here we know that the residue of the singularity originates 
in the leading term of the asymptotic expansion of the lowest order spectral function in the integrand. Therefore, 
subtracting this leading term from the spectral function itself produces a regular integral from which one can extract 
the regular series in the Laurent expansion at s — > by simple Taylor expansion. The only integral we have to do is 
then 

O s -oh + r^ rl { pi{0 " l) + ° {s) = h + > log2 - 5) + ° {s) ■ (4 - 9) 



Combining the results in Eqs. lj4~8|) and ([479]) with the fact that T(s)r(l - s) 



0(s) 



■|| in the asymptotic expansion in Eq. (|4.6I 



coefficient of the j— term in Eq. (|4.5p and hence the term 

We leave to the studious reader the pleasure of reproducing the residues of the r^prp an d 

Eq. I|4.5p . and hence the terms of O ( t^S") in Eq. (|4.6fl . using similar simple arguments. 



one easily gets the 
n^jrn singularities in 



The JZa(s) Moment. 



This moment appears in the expressions of a M 
given by the integral 



(ee/i) (eer) 



and 4 eeMAl) in Eqs. ||5T7|) . jSUSj) and (pT2T]) and it is 



© 2 ^ 2(s)= r d ^ si p2{0. 



(4.10) 



associated to the spectral function in the r.h.s. of Eq. i|2.8p . In terms of the variable S, it gives rise to the 
integral 



R 2 (s) = 2 I d8 25(1 - S 2 ) 3 - 1 h[ 1 



1 



<9 



9 



6 2 log 



1 



1 



which can be done analytically rather easily with the result 

s)(6 + 13s + 4s 2 ) T(l + s) 



2[) 9 s 2 (2 + s)(3 + s) 



r(| 



(4.11) 



(4.12) 



explicitly showing the singularities in the negative real axis. Notice that, because of the zeros provided by the 



factor 



there cannot be terms of O 



contributions to a M from two electron loops. 



/ 2 \ Pi 

( ^fi" J for p = 3/2 + n, n = 0, 1, 2, • • • in vacuum polarization 



• The fii(s) Moment. 
It corresponds to the Mellin transform (see Eq. (|3.13p ) 



^i(s) 



dio tu s 



uj ( \/4 + w 



4 + uj \i/4~+~uj + ^uj 

^2 



nW(- 



dx (1 — x) 



1-x 



no*) 



i 



(4.13) 
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This integral can also be easily done, using the expression in Eq. l|2.18p in terms of a subtracted logarithm: 



log(l - x) 



log(l - x) 



2 3 

X X 



The integral over x is then convergent at x = term by term, with the result 



(4.14) 



fii(s) = r(2-s) 



4 T(-l + 2a) 4 T(2s) 5r(l + 2s) 
3 r(l + s) + 3T(2 + s) + 9 T(3 + s) 



4T(-2 + 2s) T(-l + 2s) lT(l + 2s) 
3 = ~ 



T(s) 
4T(-2 + 2s) 

3 r( s ) 



r(l + s) 3T(3 + s) 



Hi. 



H_i 



r(-i + 2£) i r(i + 2 S ) ij 
2 r(i + s ) Hs + 3 r( 3 + s ) H2+s 



(4.15) 



where EL denotes the function 



E s = + s) + 7E . 



related to the ^-function f36(] tp(z) = ^logT(z), as follows 

4>{z) = -7e + V ( — -— - 



and, therefore, when z is an integer n, H n corresponds to the Harmonic sum 

111 1 



H„ = 1 ■ 



2 3 4 



(4.16) 



(4.17) 



(4.18) 



Notice that Oi(0) = — ^ + ^- is precisely the well known coefficient of the fourth order contribution to aju 
from the Feynman graph in Fig. 4: 



a " I J I 36 + 2 



(4.19) 




Fig. 4 Fourth-order Feynman diagram with a lowest order vacuum polarization muon-loop. 

Another property of Oi(s), which is also valid for all f2 p (s) with p > 1, is the fact that they are no longer 
singular at s = —1/2. This can be readily seen from the fact that IT^' (t^e 771 ^) m H-l^P vanishes at x = 0: 



1 — X M / a:—* \7T 



15 



x 2 + 0(x 3 



(4.20) 
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This is why Feynman graphs with electron loops and at least one muon loop insertion in the vacuum polarization 

1/21 



have no O 



\ m l) 



terms in their asymptotic contributions. However, from the asymptotic behaviour 

V4 + w — \/uj 



4 + w V y/4 + uj + ^/uj 
in Eq. (j4- 13fl . there follows that 



n w (-^ O 



a 

V 7T 



(4.21) 



fii(s) 



1 1 

'30 § + s 



+ 



(4.22) 



which induces a term of O 



of 
all s 



3/2' 



(3) 



3/2 



in the asymptotic expansion of . A priori it could also induce terms 



in a^ efl \ That this is not the case is due to the fact that as we have seen i?2(s) has zeros at 



T/2 + n for n 



explains why a term of O 



-1,-2,-3, 

/ 2\3/2 



However, as we shall later see, Rz(s) is finite at s = —3/2 and this 
does indeed appear in a{f ee/ ^ . 



The moments CIq(s) and i?2(s) fix entirely the evaluation of the contribution to a M from the sixth order Feynman 
diagram in Fig. 5, which, in our representation, is given by the integral 



c s +ioo 



/ 4»Ttj 



r(«)r(i-«) o (s) r 2 ( s ) 



(4.23) 




Fig. 5 Sixth-order Feynman diagram with two lowest order electron-loop vacuum polarization insertions. 

Since we know flo(s) and i?2(s) analytically we can evaluate the full singular series of their product and, therefore, 
by the converse mapping theorem, the asymptotic expansion of to as high a degree of accuracy as we wish. As 
an illustration, we give the results for the first few terms: 



,(«0 



1 , 2 K 



18 



mr. 



25. m' 

^r lo s 2 

54 mi 



317 

324 



27 



1/2 



45 



■log 



2 TO ' 



52 , m 
18 m 



O 



log' 



3 m l 



(4.24) 



where the contribution from the leading singularity at s — is the one in the first line, which agrees with the earlier 
calculations in refs. [ltSl. l22l|. The exact analytic evaluation of a^ e \ as well as of the full sixth-order contribution to 
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a M from electron loop insertions, including light-by-light scattering loops, can be found in the papers by Laporta and 

Remiddi [Ulll- 
On the other hand, the moments £li{s) and -ffi(s) fix entirely the evaluation of the contribution to a M from the 
mixed sixth order Feynman diagram in Fig. 6, which, in our representation, is given by the integral 



v w -(;)'<-*>Si 



ds 



( 4mi 



c s —zoo 



\ 777 



r( s )r(i- s ) n x ( s ) r^s) 



(4.25) 




Fig. 6 Sixth-order Feynman diagram with one electron-loop and one muon loop vacuum polarization insertions. 

Again, since we know f2i(s) and i?i(s) we can compute as many terms as we wish of the asymptotic expansion of 
ajf^ . We give a few terms below: 



-,( e A») 







119 



115 
"27 



9 



log 



O 



m l _ 


61 


7T 2 


ml 


162 " 




(™L 


) lof 


mi 



(4.26) 



in agreement with earlier calculations [24], [25| ■ 

Finally, we observe that once the moments fli(s) and i?2(s) are known, we have all the ingredients to evaluate the 

eighth-order contribution a£ in Eq. I|3.15p which we shall discuss in the next section. With i?2(s) known we can 

also attempt the evaluation of a^ eT ^ in Eq. I|3.17p . which we shall do in Section VII. 



The i?3(s) Moment. 



(eee) (eeeu) 

a). 



This moment appears in the expressions of et M , u, M 
given by the Mellin transform 



and a^ eeT) in Eqs. [gHH ((3T20|) and (|3~23|) and it is 



(-) Ms) = J dee-'paco. 

of the spectral function in the r.h.s. of Eq. I|2.9p . It gives rise to the (^-integral 



(4.27) 



dS 26(1 - S 



2\s-l 



3(5 



2 6 



8 1 

9 ~ 3* 

31 



log- 



(4.28) 



Only the terms proportional to log 2 i-p| require special attention. They give rise to integrals of the type 



f(a) = J\s(l-S 2 r log 2 ^ 



with a = s + n— 1 and n = 0, 1, 2, 3,4, 5 . 



(4.29) 
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Integrating by parts, we observe that /(<r) obeys a simple functional relation 

(1 + 2a)f(v) - 2a f (a 1) + h/l = Q 

from which, and the boundary condition /(0) = there follows that 

r- r(i + ff ) 



(4.30) 



r(| + a) 



■0 (1) (1 + cr) , where V (1) (z) 



dz 



V(z) 



(4.31) 



One can then do the integral in Eq. (|4.28|) . After some rearrangement so as to exhibit explicitly the singular 
structure, we get the following expression 



R 3 (s) = 



A r( a ) 

864 r(^ +s) 



Pr(s) 



- (l + s)(35 + 21s + 3s 2 ) (27tt 2 - 162 



_s(l + s)(2 + s) 
with i-V(s) the seventh degree polynomial 

P 7 (s) = 3492 - 8748s - 26575s 2 - 9214s 3 + 18395s 4 + 17018s 5 + 5120s 6 + 512s 7 . 



(4.32) 



(4.33) 



The function ^^(z) is called the Polygamma function of index one. In general, the Polygamma function of index 
n is defined as: 



ibW(z) = %b(z), with ^ 0) (z) = ib(z) for n = 1,2,3,.... 

dz n 

These functions are also related to the Hurwitz function (also called the generalized zeta function) 

oo 

((s,z) = ^T(m + z)- s 2^0,-1,-2,..., Res>l, 

as follows 



m— 



if,W(z) = (-l) n+1 n\ C(n + l,z). 



(4.34) 



(4.35) 



(4.36) 



The Polygamma functions are therefore meromorphic with poles at z = 0, — 1, —2, . . . , with multiplicities n + 1. In 
fact, the Hurwitz function and therefore the Polygamma functions, have the following Mellin representation: 



C(s,z) = 



1 



r(*) Jo 



dt t 



s-l -zt 



1 



1 - e- 



Res > 1 , Rez > . 



(4.37) 



The moments Oo(s) and i?3(s) fix entirely the evaluation of the contribution to from the eighth order Feynman 
diagram in Fig. 1(A), while the moments fii(s) and i?3(s) fix entirely the evaluation of the contribution to a M from 
the tenth order Feynman diagram in Fig. 2(B). We discuss these results in the next section. With i?3(s) known one 
can also attempt the evaluation of a|f eeT ' ) in Eq. (|3 . 23[) . which we shall do in Section VII. 

• The fia(s) Moment. 
It corresponds to the Mellin transform (see Eq. (|3.13j) ) 



duj LO 



s-l 



4 + w V ^/4 + uo + 



n("'(- t 



1 , , J x 2 

dx (1 — x) I 

i.i \1 — x 



m„ 

1-x 



(4.38) 



This integral can also be done, using the expression in Eq. (|2 . 18[) and replacing the logarithms with subtracted 
logarithms, as in Eq. (|4. 14|1 . One is then left with three types of integrals Sp\s,n) defined below, where the 
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upper index j refers to the power of the logarithm in the integrand , the lower index p to the degree of the 
subtracted polynomial and k e N : 



(4.39) 



S«( S ,n) 



Jo 



dx x Zs+n {l-x) 1 



log(l - x) 



x' 

T 



x° 

y 



x 1 

T 



x- 



60 L 1 (4 + n + s) 



30 



r(3 + n + 2s) 
r(5 + n + a) 



|2Q r(4 + n + 2.s) + i5 r(5 + n + 2.s) + i2 _r(6 + n + 2s) 



+60 



T(6 + n + s) 
r(l + n + 2s) 
r(3 + n + s) 



r(7 + n + s) 

(Hl_ s — H 2 + n +s) 



T(8 + n + s) 



(4.40) 



and 



f 1 

S^ 2) (s,n) = / da;a; 2 " +n (l-a;) 1 
Jo 



2 , / 3 11 3? 5x 13Tx 



log(l- a :)-[^ + ^ + — + — + 18Q 



- ~ 180 r(2 " S} j" 180 r(5 + n + S ) ~ 18 °r(6 + n + S ) 

165 r(5 + n + 2s) - 150 r(6 + n + 2s) - 137 r(? + n + 2s) 



r(7 + n + s) 
180 



T(8 + n + s) r(9 + n + s) 

[r(l + n + 2s) (V/(2 - s) - 2?A(2 - s)^(3 + n + s)+ 



r(3 + n + s) 

+V> 2 (3 + n + s)+ ?/> (1) (2 - s) - ?/> (1) (3 + n + s)) | . 



(4.41) 



The analytic result for ^{s) can then be expressed in terms of these three integrals as follows: 



27 

-f 5«( fll -4) + f 5«( fll -5) 



27 



+ i sf (5, 0) - ^ sf (5, -2) + sf (a, -3) + 4 S< 2) (s, -4) 



16 



^(,,-5) + ^ (,,-6). 



16 



(4.42) 



In particular, the value of ^(s) at s = fixes the contribution to the muon anomaly from the Feynman diagram 
in Fig. 7, i.e. 



943 4 2 8 . 
324 " 135^ + 3 C(3) 



(4.43) 
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Fig. 7 Sixth-order Feynman diagram with two lowest order muon-loop vacuum polarization insertions. 



Knowing r2 2 (s) and Ri(s) we have all the ingredients to evaluate ajf in Fig. 1C and knowing f2 2 (s) and Ri(s) 
we can also evaluate ajf 6 ^ in Fig. 2C. We discuss these calculations in the next section. 



The FU(s) Moment. 



This moment appears in the expression of a£f eee ' and it is given by the Mellin transform 

/(T\4 r°° 

(-) ra( S ) = j o de^-'p^o. 

of the spectral function in the r.h.s. of Eq. 1|2.10[1 . It gives rise to the <5-integral 



(4.44) 



R 4 (s) = / dd 26(1 - 6 



2\s-l 



-4tt 2 



9 3 




8 1*2 W 1 1*2 



2 6 J & l + 6 



si 1 -- 1 -** 

2 6 



log 



1-6 



(4.45) 



Here only the integrals proportional to log 3 ^rf are new with respect to the integrals which already appeared 
in the evaluation of i? 2 (s)- By integration by parts, they can be reduced to integrals proportional to log 2 £fy 
and hence to the function f(s) in Eq. (|4,31|1 . In terms of f(s) and the auxiliary function 



g(a) = 2^ 



F(a+i) 



25 



(4.46) 



we then find: 
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Ra(s) = - 



250 
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10 2 1 
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1 

54 
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(4.47) 



Knowing explicitly the function R,4,(s) will allow us to perform the evaluation of a M 
section. 



(eeee) 



which we do in the next 



• The 0,3(3) Moment. 
It corresponds to the Mellin transform (see Eq. p,13j) ) 



- n 3 («) 

7T / 



s-1 



do; w 



(1 



4 + uj \ V4 + uo + 



1 - X 



1 



(4.48) 



As in the evaluation of the integral in Eq. (|4.38p , the first step consists in replacing the logarithms and their 
powers by subtracted logarithms and subtracted power logarithms. One is then left with integrals of the type 
Sp\s,n) like those already introduced in the evaluation of 0,2(3), The new ones are: 



dx x 2s+n {l~x)^ s 
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(4.49) 
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[ dx x^il-x) 1 ' 6 [log 2 (l 
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(4.50) 



and 
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(4.51) 



The analytic result for ^(s) can then be expressed in terms of the Sp\s,n) integrals as follows: 
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In particular, the value of ^3(5) at s — fixes the contribution to the muon anomaly from the Feynman diagram 
in Fig. 8, i.e. 




Fig. 8 Eighth-order Feynman diagram with three lowest order muon-loop vacuum polarization insertions. 
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(4.53) 



Knowing fls(s) and Ri(s) we have all the ingredients to evaluate a^ t e ' w ^ in Fig. 2D. We discuss these calculations 
in the next section. 



V. EIGHTH ORDER RESULTS FROM ELECTRON AND MUON VACUUM POLARIZATION LOOPS 

We have now all the ingredients to proceed to the calculation of the eighth order contributions illustrated by the 
Feynman diagrams in Fig. 1. 

• Three Electron Loops, Fig. 1(A) [one diagram]: 
We recall the expression in Eq. (|3. 14|) 

< ee) = (^Tin J ds (tf) V-)r(i-.)no«ife(-)- (5.i) 

c s — ioo 

The converse mapping theorem relates the asymptotic behaviour of CLfj, as a function of the small mass ratio 
4TOg/m^ to the singularities of the integrand in this equation as a function of the Mellin s-complex variable. 
Using the explicit expressions for Qo(s) and i?3(s) given in Eqs. I|4.ip and (|4.28[) we can proceed to the calculation 
of the singular series expansion of the integrand in question. The relevant singularities are those in the left-hand 
side of the Mellin s-plane. They occur as multipoles at s = 0, —1, —2, —3, . . . , because of the factors T(s) and 
i?3(s); and as single poles at s — —1/2, —3/2, —5/2, . . . because of the factor f2o(s). The leading singularity at 
s — > has a quadruple pole, a triple pole, a double pole and a single pole. Their residues govern the successive 

terms of O(log 3 ^|), of C(log 2 ^f), of ©(log and of C(Cte.) in the leading contributions to a^ eee) for ^| 
large: 
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The next to leading singularity at s 
bution ofO(^): 



-1/2 is a simple pole, induced by ri ( — 1/2) , which governs the contri- 
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The next singularity at s — ► —1 in the Mellin plane is again a multipole type singularity which, by the inverse 
mapping theorem, governs the terms 
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There is no problem in evaluating as many terms as we wish. In the Appendix we give the result of the 



asymptotic contribution to a£f ee ^ up to terms of O 



( 7T ) ( mj ) 



2 \ 5/2' 



. Because of the present experimental error 



in the — - ratio, there is no need to go beyond this approximation, but we have checked that including higher 

„2\3l 



order terms up to O 



(") (m|) 



3/2' 



does not change our final numerical result in Table 1 in the Appendix. The 



terms up to O (~) agree with those obtained by Laporta [8| using a very different method. The 

agreement with the numerical determination by Kinoshita and Nio Q is quite remarkable, although our result 
is of course more precise. 

Two Electron Loops and One Muon Loop, Fig. 1(B) [three diagrams]: 
The corresponding expression is the one in Eq. (|3 . 1 5[1 
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with fii(s) and i?2(s) given in Eqs. I|4.15p and <|4. 12[) . The relevant singularities in the Mellin s-plane occur at 
s = 0, —1, —2, —3, • • • as multipoles. The singularity at s = governs the terms 
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in the asymptotic expansion. The next -to -leading singularity is at s = — 1 and it is responsible for the higher 
order terms 



(ee/i) 



(5.7) 



The analytic results in Eqs. I|5.6p and (|5.7[) agree with those given by Laporta in ref. 0. As discussed previously, 
there is no singularity at s = —1/2 because of the presence of a muon loop self-energy (see the discussion around 
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Eq. I|4.20p ). There could be a priori a singularity at s = —3/2 because it appears in the singular expansion of 
fii(s); however, it is screened by the fact that i?2(s) has a zero at s = —3/2. In the Appendix we give the result 
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for 



(ee/i) 



including terms of O 



Taking into account higher order terms brings in contributions 



which numerically are of the order of the error induced by the present value of the —£ ratio in the previous 
terms. Numerically, our result is more precise than the one obtained by Kinoshita and Nio though it agrees 
very well within their quoted error. 

• One Electron Loop and Two Muon Loops, Fig. 1(C) [three diagrams]: 
The corresponding expression is the one in Eq. (|3.16|) 
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(5.8) 



with f22(s) and Ri(s) given in Eqs. I|4.42|l and l|4.3[l . The relevant singularities in the Mellin s-plane occur at 
s = 0, —1, —2, —3, • ■ • . The singularity at s = is a double pole and governs the terms 
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(5.9) 



The next -to-leading singularity at s = —1 is a single pole which governs the next -to-leading term 
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and so on. The results in Eqs. I|5.9|) and (|5 . 10[1 agree with those of Laporta @. In the appendix we have also 



included terms up to O 



These are the terms one needs to fix to the accuracy required by 



the present knowledge of the mass ratio. It is not surprising that this analytic evaluation allows for a more 

accurate determination of this term than the numerical estimate by Kinoshita and Nio [s| , though it agrees with 
it within their given error. 



VI. TENTH ORDER RESULTS FROM ELECTRON AND MUON VACUUM POLARIZATION LOOPS 

We also have all the ingredients to proceed to the calculation of the tenth order contributions illustrated by the 
Feynman diagrams in Fig. 2. 

• Four Electron Loops, Fig. 2(A) [one diagram]: 
We recall the expression in Eq. I|3.19p 

c s +ioo 

a ( eeee) = ^ 5 J ds(^pj T{s)T{l-s)n Q {s)R A {s), (6.1) 

c s — ioo 

with Qq(s) and Ri(s) given in Eqs. (|4.ip and (|4.47p . The converse mapping theorem relates the asymptotic 
behaviour of flu ctS 3b function of the small mass ratio 4r7ig/m^, to the singular series expansion of the 
integrand in this equation. The relevant singularities are those in the left-hand side of the Mellin s-plane. 
They occur as multipoles at s = 0, —1, —2, —3, . . . , because of the factors T(s) and R^s); and as single poles 
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at s = -1/2,-3/2,-5/2, 
contributions: 



because of the factor Qq(s). The multipoles at s — > govern the asymptotic 



(eeee) 




162 



log 



4 m l 



/8609 
V4374 
64613 
" 26244 + 729 7 



50 , 

7T 

243 
317 „ 



25 3 ml 
" 243 l0g ^ 

8 \ m? 
27 / m; 



135 



100 
~87 



C(3) 



317 

486 



81 



log 



2 



(6-2) 



There is only a simple pole at s = —1/2. It governs the next-to-leading asymptotic contribution 
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The contributions in Eqs. (|6.2|1 and l|6.3p agree with those given by Laporta in ref. [0]. We have calculated 

m 2 

further contributions, up to the required accuracy governed by the present knowledge of the ratio. They are 

given in the Appendix. Numerically, although our result is much more precise than the one given by Kinoshita 
and Nio 0, it agrees very well with it within their quoted errors. 

• Three Electron Loops and One Muon Loop, Fig. 2(B) [four diagrams]: 
We recall the expression in Eq. (|3.20p 
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with fii(s) and i?3(s) given in Eqs. I|4.15p and (|4.28p . The relevant singularities here occur as multipoles at 
s = 0, — 1, —2, —3, . . . , because of the factors T(s) and i?3(s); and as single poles at s = —3/2, —5/2, . . . because 
of the factor f2i(s). The multipoles at s — > govern the asymptotic contributions: 
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This expression agrees with the one given by Laporta [§]. We have also calculated the corresponding terms 
f / 2 \ 5/21 

up to 



( m l ) 



as required by the wanted accuracy. Again, although numerically our result is much more 
precise than the one given by Kinoshita and Nio @|, it agrees very well with it within their quoted errors. 

• Two Electron Loops and Two Muon Loops, Fig. 2(C) [six diagrams]: 
We recall the expression in Eq. I|3.2ip 
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with (s) and i?2(s) in Eqs. <|4.42|) and (|4.12p . Here, the multipoles at s = generate the terms 
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which agree with the expression found in Laporta [9J], We have also calculated the contributions up to 

„2\3 







3 % 



, which can be found in the Appendix. Our numerical result agrees with the one given 
by Kinoshita and Nio [H] within their quoted errors. 

One Electron Loop and Three Muon Loops, Fig 2.(D) [four diagrams]: 
We recall the expression in Eq. (|3.22[) 



,2' 



with ^(s) and Ri(s) in Eqs. l|4.42p and (|4. 12[) . Here, there is a double pole at s = which generates the terms 
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in agreement wit the expression found in Laporta 0|. We have also calculated the contributions up to 

„2\3 



o 



, which are given in the Appendix. Our numerical result also agrees very well with the 



one given by Kinoshita and Nio [j| within their quoted errors. 



VII. CONTRIBUTIONS FROM ELECTRON LOOPS AND ONE TAU LOOP 



We shall finally discuss the calculation of vacuum polarization contributions involving electron loops and a tau loop. 
The corresponding expressions are the ones in Eqs. I|3.17p and (|3-23|1 . As already mentioned, the problem here is the 
non factorization of the dependence in the the two Mellin variables s and t because of the presence in the integrand 
of the function 9(s, b) defined in Eq. I|3.18p . 



A. Two Electron Loops and One Tau Loop 



The corresponding expression in this case is the one given in Eq. p,17p . Using the result obtained for i?2(s) in 
Eq. I|4.12p and the expression for 0(s, t) in Eq. I|3.18p . we have: 
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The evaluation of the asymptotic behaviour of this type of integral calls for a more sophisticated material than 
the inverse mapping theorem applied to the previous cases where there was only one ratio of masses. Indeed, the 
generalization of the inverse mapping theorem to this case is a typical problem of calculus of residues in C 2 . Let us 
then adopt the standard notation of multidimensional complex analysis [33] and denote by o/ eer ) the following 2-form: 



JeeT} = /4m|V s / rr^ \ _< T(t) T(l-t) r 2 (2 - t) ^ 




V m 2 J \m 2 T j t r(4 - 2t) 
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s 3 (2 + s)(3 + s) r(s + |) 

r(i + 2s - 2t) r(2 - s + 1) , , , 

r (3 + s-t) ~ dsAdt - ^ 

Recall that the Mellin-Barnes representation in Eq. (|7,ip is valid for all s and t such that Re(s) G]0, 1[ and Re(t) £ 
] — 1, 0[; these two conditions resulting in the grey fundamental square in the plane [Re(s), Re(f)], as illustrated in 
Fig. 9. In the case of two Mellin variables, this fundamental square, which in general may become a fundamental 

polyhedra, generalizes the concept of the fundamental strip in the case of one variable. Since both and are 

small, the asymptotic behaviour we are looking for in Eq. (|7.ip . will be governed by the residues associated to the 
singularities of the cone jRe(s) < 0, Re(i) < —1}, denoted by II in Fig. 9. Formally, the solution to our 

problem is then given by the sum [38] 



a< r ) = (%y E R«wo)" (eeT) . (7-3) 

(s ,to)en 

From the expression in Eq. I|7.2p . one can easily see that the singularities in the plane [Re(s), Re(t)] associated to 
w' eeT ', are defined by a set of straight lines. For example, T(2 — s + t) induces a family of singular lines parameterized 
by the affine equation: Re(i) = Re(s) — (2 + n), n > 0. Each of the singular lines is called a divisor. As discussed in 
ref. [2§| , it is sufficient for our purposes to consider the singular set of points defined by the intersections of all the 
divisors in the appropriate cone. In our case, a subset of the nearest singular points to the origin of the cone II, is 
illustrated by the red dots in Fig. 9. 
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Fig. 9 Plot of the divisors and the singular set of points of or 667 ^ in Eq. The small red dots are singularities of the first 

kind; the larger red dots are singularities of the second kind. Notice that the singularities on the line s — — 2 , k > 1, are 
screened by the presence of the factor T(s + |) in the denominator of ui eeT \ 

There are in fact two kinds of singularities to consider: 

1. Singularities of the first kind, where only vertical and horizontal divisors (two or more) cross each other. 
In order to obtain the residue associated to such a singular point (so, to), one proceeds as follows: 
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• First perform the change of variables s <-> so + s , t «-> to + f so as to bring the singularity to the origin 
(0,0). 

• Use the functional relation T(z+ 1) = zT(z) for each singular Gamma function until reaching an expression 
of the type 

w H = !M dsAdt] ( 7 .4) 

s n t m V > 

explicitly showing the full singular behaviour at the origin, with h(s,t) an analytic function at (0,0) and 
n and m positive integers. 

• Then, use the Cauchy formula 



~ , leer) 1 d n+m ~ 2 h( S ,t) 
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2. Singularities of the second kind, where two or more divisors cross each other, one or more of them being 
oblique lines. 

This is the case, for instance, at the point (0, —2) in Fig. 9. In fact, in our case, the singularities of the second 
kind at a point (so,io) are induced either by the oblique divisors generated by T(l + 2s + 2t) or the oblique 
divisors generated by T(2 — s + t). In order to obtain the residue associated to a singularity of the second kind 
at a generic point (so>*o), one proceeds as follows: 

• Again, first perform the change of variables s «-* So + s, t <-> to + 1 so as to bring the problem back to 
a singularity at the origin (0,0); and then apply the functional relation T(z + 1) = zT(z) to the singular 
Gamma functions in question, until able to write o/ eer ) in a way which explicitly shows the full singular 
behaviour in question: 

Jeer) = ^M) ds A dt , (7.6) 

with h(s, t) an analytic function at (0, 0) and n and m positive integers. The factor t — s in the denominator 
is specific to the class of oblique divisors in our case which are all parallel lines to s = t with multiplicity 
one. 

• Apply the Transformation Law for Residues [26| , to the form 

Jeer) m 

Eq. I|7.6p so as to fully decouple the 
s and t dependence in the denominator. This, of course, requires some explanation which we next provide. 

To be precise, let us discuss the case of the singular point (0, —2) in Fig. 9. There is a very important theorem in 
multidimensional complex analysis, known as the Transformation Law for Residues [2(| which, in the case of C 2 , 



states that if U is an open set containing (0,0); if f 



and g 



9x(s,t) 

32 (M) 



are two analytic mappings 



from U to C 2 such that f 1 (0, 0) = g 1 (0, 0) = (0, 0); and if there exists an analytic two- by-two matrix A such 
that g = Ai , then 

h(s,t) , , ^ h(s.t) detA(s.t) , . x 

/i(M) h{s,t) 9i{s,t) 52 (M) 

The application of this theorem to our case goes as follows: 

• First note that with the change of variables s — * s — t, t — > t, 

h(s,t) , , h(-s + t,t) , . , 

Res (0 , 0) ds A dt = Res (0 , 0) (is A dt . (7.8) 

We do this because in this form and as shown in the next step, we can then find easily a matrix A to apply 
the Transformation Law for Residues. 
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Indeed, with f the vector: 



st 

(-s + t) 3 

we can find a matrix A which does the transformation: 

-t 2 + 3si-3s 2 a 
s 2 - 3st + 3t 2 t 



'-s 4 ' 
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=g 
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and 



detA = -(s 3 + t 3 ) . 

• By virtue of the theorem quoted above, we can then assert that 
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(7.12) 



where in going from the second line to the third, we have applied the Cauchy formula in Eq. (|7.5|l . 
Finally, from Eq. (|7.8p , and applying the formula for chain-derivation in Eq. (|T. 12|) , we get the result 
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We shall now proceed to the calculation of a^ eer ^ in Eq. (|7.3p . following the procedure which we have just outlined. 

• Singularities on the line t = — 1. 

- The leading singularity is at the point (so, to) = (0,-1), and it is a singularity of the first kind. After 
translation to the origin and the reduction to an explicit singular form, we can write o/ eeT ' in the following 
way: 



(eer) = V-i)(a,t) 
s 3 t 



ds A dt , 



(7.14) 



with 
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Then, using Eq. 17. 5^ we get 
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— The next singularity on the line t = — 1 is at the point (so,to) = (—1,-1), and it is also a singularity of 
the first kind. In a similar way to the singularity in (0, —1), one can easily find that the residue associated 
to this singularity is given by the expression 



which gives a small contribution because of the suppression factor 




• Singularities on the line t = — 2 
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— log- 
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The singularity nearest to the origin on the line i = — 2 is at (s , t ) — (0, —2) and it is of the second kind. 
After translation to the origin and the reduction to an explicit singular form, we can write o/ eer ) in the 
following way: 
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where 
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Then, using the result in Eq. I7.13l we get 
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which for convenience we have ordered in decreasing powers of ( — f ) log™ — f with n = 3, 2, 1, 0. 

The other singularities on the line t = — 2 (for s = —1, —2 and further) turn out to give contributions which are 
too small to be of physical relevance. 

• Other Singularities 

In fact there is only one more singularity which can give rise to a term of the order of the present experimental 
error limitations. This is the singularity at (0,-3). Its residue can be computed exactly in the same way 

( 2 \ 3 ^2 
\ g n ^ with 

n = 3,2,l,0: 
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The total contribution from Eqs. l|7.16p . I|7.17p . (|T.20|) and (|7.2ip to a^ eT \ ordered in a decreasing order of 
magnitude, is given in Eq. (|3. lTfl in the Appendix. 



B. Three Electron Loops and One Tau Loop 



The corresponding expression in this case is the one given in Eq. (|3.23p . Using the result obtained for R%{s) in 
Eq. I|4.28p and the expression for 0(s, t) in Eq. (|3 . 18[) . we can write explicitly the two form associated to this integral 
as follows: 
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r(3 



t) 



ds A dt . 



(7.22) 



where P>r(s) is the polynomial given in Eq. (|4.33|1 . When comparing the two forms J^ eT ^ in Eq. I|7.2p and w 
in Eq. (|T.22|1 . we can see that they only differ in the form of the factorized s-dependence. Therefore, except for the 
multiplicity of the vertical divisors, the plot of the singular points associated to uj^^ which we show in Fig. 10 is 
pretty much the same as the one for oj^ eeer ' in Fig. 9. It is then not surprising that the calculation of the asymptotic 
behaviour of the integral in Eq. (|3.23|) turns out to be very similar to the one discussed in the previous subsection. 
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Fig. 10 Plot of some of the divisors and singular points of w/f 6 with the same conventions as in Fig. 9. Only the 
multiplicity of the vertical divisors differs from the former case. 
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Residues at the singular points in Fig. 10 are computed following the strategy discussed in the previous subsection, 
where one uses the functional relations: T(z + 1) = zT(z) and ip^(l + z) = ip^'^z) — \. 

For the singularities of the first kind one can then use Eq. I|7.5p . The new type of singularities of the second kind are 
reducible to a form of the type 



(eeer) 



h(M) 



s i t(-s + t) 



ds A dt . 



(7.23) 



with h some analytic function at the origin. In this case, the appropriate A-matrix to implement the Transformation 
Law for Residues reads: 



-t 3 + 4si 2 - 6s 2 t + As 3 s 
-s 3 + As 2 t - 6st 2 + At 3 t 



and proceeding as explained in the previous subsection one finally obtains: 



h(s,i) , , 1 
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rd 4 h(M) , ^(M) | fi d 4 h(M) 
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dsdt 3 
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d 3 hs,t) 
ds 3 dt 



(7.24) 
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(0,0) 



Let us now discuss the explicit calculation of a£f eer ^ . 



• Singularities on the line t = — 1. 

- As in the case of a^ eT \ the leading singularity is the one at (so>*o) = (0, —1) and it is a singularity of the 
first kind. Performing the change of variables t «-> — 1 + t we write: 

(eeer) _ h (0,-l)M) 



or 



sH 



■ ds A dt , 



where 
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From this expression we obtain 
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which provides the leading contribution to a)f e 

The next singularity is at (so, to) = (—1, —1). Also a singularity of the first kind, whose residue we find to 
be 



Res ( _i,_i) w (eeeT) = 
and gives a contribution suppressed by a factor 



4 , 2 to 2 20 , m 2 634 8tt 2 
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Singularities on the line t 



-2. 



— The singularity at (sq, to) = (0, —2) is of the second kind. After the change of variables t <-> — 2 + t, and 
using standard functional relations, one can write o/ eeeT ) in the following way 



„(^ ) = V- a) (*,t) daAdt 



s 4 i(-s + 1) 



(7.30) 



with 
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According to Eq. l|7.25p . we have: 
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^ log™ 3 with n = 4,3,2,l,0. 

• Singularities on the line t = -3. 



The singularity at (0, —3) is also of the second kind, and its residue can be computed in the same way as 
the one for the singularity at (0, —2). One finds: 
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Res (0 ,_ 3) ^ = -f 
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446596444350000 1913625 42525 893025 

( 2 \ 3 ^2 
— f ) log" — f with n = 4, 3, 2, 1, 0. 

The other singularities have residues which give contributions smaller than the error induced by the leading 
term. The total contribution from Eqs. I|7.28p . (|7.29[) . i|7.32p and <|7.33[1 to aj?f e , ordered in a decreasing order 
of magnitude, is given in Eq. (|B.9[) in the Appendix. Notice that, numerically, the contribution with one tau 

loop and three electron loops: a^ eeT \ is of the same size as the contribution from three muon loops and one 
electron loop: . 

APPENDIX 

The purpose of this appendix is to collect systematically the various results discussed in the previous sections. This 
may be useful to readers who are only interested in the final analytic expressions and the numerical results. The 
values that we have used for the lepton masses are the ones in the 2006 PDG booklet (30l ] : 

m e = 0.510 998 92(04) MeV , m„ = 105.658 369(9) MeV and m T = 1776.99(29) MeV . (0.34) 



A. Eighth Order Results 

These are the results corresponding to the Feynman diagrams in Fig. 1 with the combinatoric factors included. The 
numbers in italics are the results from the numerical evaluation of Kinoshita and Nio in ref. 0|. 
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B. Tenth Order Results 



These are the results corresponding to the Feynman diagrams in Fig. 2 with the combinatoric factors included. 



Our numerical results for 



(eeee) (eee/i) (ee/i/^) 



and o£ wx/ ^ agree, within errors, with those of Laporta || which, 
however, he obtained using an older determination of the mass ratio ^ = 206.768 262(30). The number in italics 
are the results from the numerical evaluation of Kinoshita and Nio in ref. [6j]. 
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